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Zero-field muon spin relaxation (ZF-/^SR) data for dilute spin magnetic systems have been widely 
interpreted with what is called a Kubo-Toyabe form based on a Lorentzian distribution of local field 
components. We derive here the proper magnetic field magnitude distribution using independent 
and uncorrelated component distributions. Our result is then compared to the previously accepted 
formula for ZF-/iSR. We discuss the origins of the magnetic field component and magnitude distri- 
butions. Further we found that after rescaling the magnetic field, the differences that are amenable 
to experimental examination are quite small, although the interpretations behind them are quite 
different. 
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Zero- field Muon Spin Relaxation (ZF-/iSR)[l] has been 
used long and widely to study the local magnetic field dis- 
tribution in the rare-earth metallic alloys, spin glasses, 
heavy fermion systems, superconducting compounds and 
other magnetic systems. [2] R has been believed[3] that 
in the case of the dilute limit of sparse magnetic mo- 
ments, the local field component distribution is Cauchy 
or Lorentzian-like: 
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where Bi = B^, By or B^, and a is the half- width at 
half- maximum (HWHM). 

On the other hand, from studies of systems with ran- 
domly distributed dilute magnetic impurities it was con- 
cluded that the internal field magnitude distribution is 
expected to be: [4-7] 
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Consequently the ZF-/iSR relaxation function takes the 
Lorentzian Kubo-Toyabe form 
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where A = 7;^r, and 7^ is the gyromagnetic ratio of the 
muon 27r(13.554 kHz/G). 

In the following calculation, we will prove that the 
magnetic field magnitude distribution corresponding to 
the component distribution (1) is: 
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R the components of the magnetic field are independent 
and uncorrelated with each other then the three dimen- 
sional probability distribution p(B) = g{Bx)g{By)g{B z) ■ 
Note that this probability distribution is not isotropic. 
In fact the three dimensional distribution as a product of 
individual component distributions is only isotropic for 
the case of gaussian individual component distributions. 

The magnitude distribution Pc(|B|) may be written as 
an integral over the p(B): 
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Inserting Eq. (1) into the expression of Pc(|B|) above. 
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In the equation above, we integrate (/? first. By using 
the residue theorem 
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Putting this result back into Eq. (7), we get 
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By using the identity Eq. (6), we finaUy get the mag- 
netic field magnitude distribution Eq. (4). A study of 
this field distribution discloses two important features: 

a) At low field, the field magnitude distribution, Eq. 
(9), is asymptotic to ^-t-, i.e. proportional to |B| . 
This is the same asymptotic behavior shown by Eq. (2), 
a result expected for a variety of models. [4] However, 
it should be noted that it is difficult to obtain a spe- 
cific heat which is linear in T at low temperatures as 
observed experimentally [8] with a field magnitude dis- 
tribution which is proportional to |B| while using the 
molecular-field model with the spins of the impurities 
quantized along the local field and no overall preferred 
spin direction. [4] On the other hand, there exist other 
physical models [9, 10] which would give at low field a con- 
stant field-magnitude asymptote. These models would 
yield a linear low temperature T dependence for the spe- 
cific heat. 

b) It can be determined that the Wronskian of Eq. (4) 
and Eq. (2) (even the rescaled Eq. (2)) is nonzero which 
implies that they are indeed independent i.e. Pc(|B|) ^ 
C • P(£' • |B|) where C and D are arbitrary constants. 
Nevertheless, it is interesting to compare these two dis- 
tributions numerically. If one starts from Eq. (2) and 
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1 and 



h{v) and f{u) have maxima at vq ~ ^ 

IR I 

uq = - — - ~ 1.476 respectively. In order to compare 
these two distributions we let them reach maximum at 
the same magnetic field |B|. Then we obtain the follow- 
ing relations: 
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where j{u) is just a rescaled distribution of h{v) for u 
1.476V. 



(where t ~ cost 
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FIG. 1. Comparison of the probability distributions of the 
magnitude of the magnetic field between f{u) and j{u). 

A plot showing the similarity and difference between 
f(u) and j(u) is shown in Fig. 1. 

It is interesting to study Eq. (2) further. Suppose we 
divide it by 47r |B| to change it to a probability density 
distribution (we assume here the field is isotropic) and 
do a double integral in the Cartesian coordinate system 
as follows: 
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If we take F = a, the above equation is exactly Eq. 
(1). At first glance, this is quite puzzling since we started 
from the same distribution and obtained a totally differ- 
ent magnitude distribution: Eq. (4). A careful study of 
all the derivation procedures above, leads us to conclude 
that: there is not a one to one mapping or correspondence 
between the component distribution and magnitude dis- 
tribution; to some extent, the two seemingly paradoxical 
results are due to different causalities; the three com- 
ponent distributions behind the magnitude distribution 
Eq. (4) are assumed to be independent and uncorre- 
lated, on the other hand, the component distribution 
Eq. (14) arises from the magnitude distribution assum- 
ing isotropy. 

It can be seen from figure 1 that although the form of 
our new result is analytically quite different from the pre- 
vious one which is widely used in ZF-/iSR spectroscopy. 



numerically they are amazingly close to each other. As 
far as the relaxation function is concerned, an experimen- 
tal determination of which field magnitude distribution 
is correct should be essentially impossible. 

A commonly- used procedure [4] which yields the 
Lorcntzian field magnitude distribution of Eq. (2) is to 
assume a RKKY(Ruderman-Kittcl-Kasuya-Yosida) in- 
teraction among the dilutely distributed spins and to use 
the MRF(Mean Random Field) approximation. Further 
this produces the distribution of the vector field B (not 
just the individual component distributions). Even for an 
isotropic vector field distribution, it rarely occurs that 
this distribution factors into 3 uncorrelated field com- 
ponent distributions. Our approach to this problem is 
logically straightforward and independent of any phys- 
ical models, the only assumption here is the use of an 
appropriate field component distribution. 
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